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基于 分 层 介质 模型 的 HF 雷达 高 度 计 
火星 回 波 仿真 研究 ” 
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(中 国 科学 院 国家 空间 科学 中 心 微波 遥感 技术 重点 实验 室 ”北京 100190) 


摘 要 火星 次 表层 中 的 水 探测 是 当前 深 空 探测 领域 中 的 热门 问题 之 一 . 高 频 (HF) 雷达 高 度 计 具有 较 强 的 穿 透 
能 力 , 且 能 同时 实现 距离 和 功率 测量 , 是 火星 次 表层 探测 的 重要 手段 . 介绍 了 高 频 雷 达 高 度 计 系统 的 原理 和 设计 ， 
通过 分 析 高 度 计 电磁 脉冲 与 多 层 光 滑 介 质 之 间 的 相互 作用 , 得 出 其 回 波 功率 随时 间 变 化 的 模型 , 并 考虑 了 表面 粗糙 
度 对 表面 回 波 功率 的 影响 ， 选 用 典型 的 火星 分 层 介质 模型 ， 对 其 介 电 常数 特性 和 高 度 计 回 波 波 形 进行 了 仿真 . 仿 
真 结果 表明 , 采用 高 频 雷 达 高 度 计 系统 可 实现 对 火星 次 表层 介 电 特 性 垂直 廓 线 的 反 演 , 对 火星 中 水 的 识别 具有 重要 
作用 . 

KE) 高 频 畦 达 高 度 计 , KE, 分 层 介 质 模 型 , 介 电 常数 , 回 波 功率 
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Abstract The detection of water in Mars is currently one of the most striking issues in deep space 
detection. High Frequency (HF) radar altimeter has good penetrating ability and can measure the 
range and power simultaneously. In this paper, the principles and design of the altimeter system were 
presented, the interaction between the altimeter electromagnetic pulses and multi-layer smooth media 
was deduced, and the model of the echo power with respect to the time delay was extracted. Besides, 
the effects of Mars surface roughness were also considered. Finally, two typical Mars multi-layer 
models were adopted, and the dielectric constants and the altimeter echo waveforms were simulated. 
Results of the simulation show that the vertical dielectric constant profile of Mars subsurface could 


be retrieved by HF radar altimetry, and the water in Mars could be identified in certain cases. 


* 中 国 科学 院 国家 空间 科学 中 心 重点 培育 基金 项 目 资助 (Y32114FB9S) 
2014-08-11 收 到 原稿 , 2015-06-30 收 到 修 定稿 


E-mail: xuxiyu@mirslab.cn 


RES 等 : 基于 分 层 介质 模型 的 HF 雷达 高 度 计 火 星 回 波 仿真 研究 747 


Key words 


Echo power 


0 引言 


火星 中 水 的 探测 是 当前 深 空 探测 领域 中 的 热门 
问题 之 一 . 大 量 研究 表明 , 火星 表层 由 风 沉 积 层 、 硬 
沉积 层 等 构成 的 侍 壤 (regolith) 覆盖 , 水 分 子 以 固态 
或 液态 的 形式 存在 于 次 表层 区 域 1-2. 高 频 (HF) 频 
Bz (3~30 MHz) 电磁 波 具 有 较 强 的 穿 透 能 力 , 特别 适 
用 于 对 火星 次 表层 的 探测 3. 本 文 提出 HF 雷达 高 
度 计 可 实现 对 火星 表层 和 次 表层 介 电 特 性 的 探测 . 


1 HF 雷达 高 度 计 原理 和 设计 


用 于 火星 表层 /次 表层 探测 的 HF 雷达 高 度 计 是 
一 台 天 底 点 指向 雷达 , 通过 发 射 脉冲 及 测量 回 波 传输 
时 间 进 行距 离 测量 , 并 通过 调整 接收 机 的 自动 增益 控 
fi] (Automatic Gain Control, AGC) 进行 回 波 功率 
测量 I, 可 实现 对 火星 表层 尘 壤 厚 度 和 次 表层 的 
分 层 介 电 常数 廊 线 等 参数 的 反 演 , 进而 识别 火星 地 下 
AK, 并 绘制 次 表层 沿 轨 向 散射 回 波 的 剖面 分 布 图 . 

根据 科学 目标 , 参考 欧洲 MARSIS (Mars Ad- 
vanced Radar for Subsurface and Ionospheric Soun- 
ding) 系统 [6-7] 和 美国 SHARAD (SHAllow subsur- 
face sounding RADar) 系统 B- 的 载荷 设计 经 验 , 本 
文 设计 了 HF 雷达 高 度 计 , 其 系统 参数 列 于 表 1. 

为 了 在 穿 透 深 度 与 垂直 方向 分 辨 率 之 间 实 现 性 
能 折 中 ，HF 高 度 计 采用 了 20MHz 的 中 心 工 作 频 
率 . 为 提高 其 距离 分 辨 能 力 ， 高 度 计 采用 线性 调频 
脉冲 (Chirp 信号) 作为 发 射 脉冲 ,时 宽 100 ps, 带 
宽 10MHz, 压缩 比 1000. 接收 机 采用 全 去 斜 技术 实 
现 脉冲 压缩 外， 

根据 最 新 的 高 度 计 体制 , 高 度 计 采 用 了 传统 和 合 
成 孔径 两 种 工作 模式 , 该 体制 已 成 功用 于 Croysat-2 
等 地 球 遥 感 测 高 卫星 1. 合成 孔径 工作 模式 可 显著 
提高 高 度 计 的 沿 轨 向 分 辩 率 及 有 效 信 噪 比 . 两 种 工作 
模式 之 间 根 据 需 要 可 自由 切换 .根据 系统 多 普 勒 带 
宽 的 需要 , 高 度 计 设计 了 600 Hz 的 脉冲 重复 频率 . 

由 于 高 度 计 工作 波长 较 长 , 采用 半 波 侦 极 子 天 


High frequency radar altimeter, Mars, Multi-layer media model, Dielectric constant, 


RER M. 与 普通 偶 极 子 天 线 相 比 , 半 波 侦 极 子 天 
线 的 阻抗 特性 很 容易 实现 谐振 ,在 方向 图 上 无 副 辩 . 
半 波 偶 极 子 的 天 线 长 度 为 波长 的 1/2, 考虑 到 雷达 
采用 的 线性 调频 信号 工作 频率 最 低 为 15 MHz (对 应 
于 20m 波长 ), 天 线 孔 径 不 能 低 于 10m. 天 线 孔 径 较 
大 , 因此 设计 为 可 伸展 结构 . 


2 HF 雷达 高 度 计 回 波 模 型 


根据 研究 , 火星 表面 上 层 更 多 的 是 由 于 火山 、 冲 
积 、 风 、 冰 川 等 过 程 形成 的 地 层 , 在 大 尺度 上 , 该 层 
可 被 解释 为 具有 指数 衰减 特性 的 多 和 孔 风 化 土壤 , z 值 
深 处 的 孔 性 为 PA 


®(z) = 有 (0) exp(—z/K). (1) 


其 中 , 8(0) 为 表面 的 孔 性 (通常 在 0.2~0.5 内 ), 即 混 
合 物 中 空气 、 水 、 冰 等 的 比例 ( 按 体积 算 ), K 为 衰减 
速度 (采用 距离 量 纲 ). 

火星 近 表 层 土壤 的 孔 中 可 能 含有 大 量 冰 态 水 , 在 
较 深 层 则 可 能 表现 为 液态 . 图 1 给 出 了 火星 地 层 模 
型 示意 , 包括 全 球 水 和 冰 的 分 布 . 这 只 是 一 个 典型 分 


表 1 HF 频段 雷达 高 度 计 系统 参数 表 
‘Table 1 System parameters of HF-band 


radar altimeter 


工作 频率 /MHz 20 

飞行 高 度 /km 400 

信号 形式 线性 调频 信号 
AY FE / us 100 

带宽 /MHz 10 

发 射 功率 /W 100 (传统 ); 20 (SAR) 
脉冲 重复 频率 /Hz 600 

天 线 体制 半 波 偶 极 子 
天 线 口径 /mm 10 

天 线 波 东 宽度 /(*) 77 

天 线 增益 /dB 2.14 
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二 


km 
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图 1 火星 次 表层 模型 


Fig. 1 Mars subsurface model 


radar altimeter 


| so, Ho ko 
atmosphere layer | 


图 2 火星 次 表层 分 层 


Fig.2 Mars subsurface multi-layer sketch 


布 示例 , 以 下 将 通过 具体 示例 对 火星 次 表层 结构 进行 
仿真 . 
2.1 平坦 界面 次 表层 模型 

如 前 所 述 , 火星 表面 以 下 基本 都 是 分 层 分 布 的 . 
因此 首先 分 析 平 坦 界面 的 分 层 模型. 

考查 如 图 2 所 示 的 简化 分 层 模 型 . 这 里 将 其 分 
为 了 三 层 . (1) 大 气 层 , 高 度 为 Ho; (2) 沉积 层 , 高 度 
为 H; (3) ZRAZ. 设 三 层 的 电气 特性 如 下 - 

(1) 大 气 层 介 电 常数 为 so, 损耗 正切 为 tando (24 
为 0), 磁 导 率 为 po. 

(2) 沉积 层 介 电 常 数 为 sl, 损耗 正切 为 tan61, 磁 
导 率 为 Hi- 
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en; 


、 经 


(3) 玄武 岩层 介 电 常数 为 sz, 损耗 正切 为 tan62， 
HERA H2 

除 少数 铁 磁性 介质 外 , 一 般 介 质 的 磁 导 率 均 接近 
于 uo; 通常 假设 ji © po; pa © uo. 

高 度 计 发 射 的 信号 穿 透 大气 层 达到 火星 表面 , 共 
产生 两 个 回 波 (二 次 以 上 反射 的 回 波 不 考虑 , 因为 其 
能 其 足够 低 ). 一 是 火星 表面 的 直接 反射 回 波 ; 二 是 透 
射 入 火星 表面 沉积 层 的 电磁 波 经 过 衰减 和 玄武 岩 界 
面 的 反射 , 在 沉积 层 和 大 气 层 传输 后 , 由 高 度 计 接收 
的 回 波 . 

这 两 部 分 回 波 可 利用 延 时 区 分 . 两 个 回 波 之 间 的 
延 时 差 为 


可 得 
t 
en A 


这 样 即 可 根据 延迟 计算 高 度 . 由 式 (3) 计算 高 度 的 前 
提 是 需要 知道 沉积 层 的 介 电 常数 e. 该 值 可 由 表面 
回 波 功率 得 到 . 

镜面 的 后 向 散射 系数 完全 相干 , 其 雷达 方程 表现 
为 Friss 传输 方程 , 即 


_ BXG2ITE 


tr 167 万 2 ` (4) 
AF, 了 为 两 个 介质 之 间 界 面 的 反射 系数 , 有 
_ yE- yE 
r= VEE SE (5) 


AH, sl 和 &2 分 别 为 两 层 介质 的 介 电 常 数 . 
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由 式 (4) 直接 可 得 大 气 层 和 沉积 层 界面 的 反射 
回 波 功率 为 ee 
4 0 
Po = = ieee (6) 
RP, o 为 大 气 层 和 沉积 层 界面 的 反射 系数 .沉积 
层 和 玄武 岩层 界面 的 反射 回 波 功率 为 
_ PAXG2(1 = |D] 


Pa = 16r(Ho + Hy)2L? (7) 


回 波 包 括 一 次 反射 、 两 次 透射 和 两 次 衰减 , (1 一 
[Tol?) 为 电磁 波 在 大 气 层 和 沉积 层 界 面 处 的 透射 率 ， 
T 为 沉积 层 与 玄武 岩 界面 处 的 反射 系数 , L 为 电磁 
波 在 沉积 层 中 的 衰减 , 有 


m 
nfe Hi| 5 
ps fea cl _ nftanðHiy uE, _ 


C C 


uf tanóH Vi ep) = 


C 


8.697 f tent ew es (dB) ~ 


(9.11 x 1078) f tané Hi Ve’ (dB). (8) 


其 中 , c 为 真空 中 的 光速 ,cf 和 分别 为 相对 介 电 党 
数 和 磁 导 率 , u 设 为 1 式 (8) 中 ， 

1Np = 20 lg edB = 8.69 dB. 

因此 ,根据 两 个 界面 回 波 的 功率 和 相对 时 延 可 以 
确定 两 层 介质 的 介 电 特性 . 对 于 多 层 平 坦 介质 , 可 以 
按照 同样 的 方法 建立 各 层 反射 回 波 模型 . 
2.2 “粗糙 表面 的 几何 光学 模型 

2.1 节 中 假定 火星 表面 是 平坦 的 , 但 这 通常 不 符 
合 实际 情况 . 以 下 将 建立 火星 粗糙 表面 模型 . 火星 的 
表层 呈现 大 尺度 和 小 尺度 六 糖度 . 由 于 本 文 研究 的 是 
火星 轨道 器 的 吕 感 测量 ,大 尺度 模型 已 能 满足 需求 
对 于 着 陆 器 或 巡视 器 的 次 表层 雷达 , 则 需要 考虑 小 尺 
度 模型 

火星 粗粮 表面 散射 元 的 高 度 概率 密度 满足 零 均 
值 高 斯 分 布 bal, 即 

2 

an (9) 


ep (- a7 


f(h) = 


~ 2Xon 


AP, on 为 本 地 高 程 分 布 的 均 方 根 值 ， 本 地 均 方 根 
斜率 可 由 本 地 均 方 根 波 高 和 相关 长 度 ! 导出 , 即 


202 
o? = oR|Cz(0)| = =. (10) 


大 尺度 粗糙 表面 的 散射 场 主要 成 分 是 镜面 E 
镜面 ) 反射 , 在 0.2~10km 的 尺度 上 ,os 一 般 不 超 
过 0.02rad (1°). 

根据 经 典 基 尔 霍 夫 近似 下 的 准 镜面 几何 光学 模 
型 , 回 波 的 散射 系数 为 


0 [To =) 
6) = ; 
a 202 cost 0 E ( 202 


(11) 


其 中 , 0 为 电磁 波 入 射 角 . 设 偏 离 天 底 方向 的 斜 距 
为 R, 则 由 几何 关系 , 有 


cos = 2, tang = yi He (12) 
将 式 (12) 代入 式 (11) 可 得 
_ \To|?R4 R? — H2 
PS 202 H4 exp ( — HoF ) (a) 
显然 , 传输 距离 R 是 时 间 t 的 函数 , 即 
ct 
R = 9? 
ha t (14) 


VA to 作为 时 间 轴 的 零点 , Be tp = t—to (tp < to), 
则 有 


R- H= (ye -Bx a (15) 
代入 式 (11) 并 化 简 可 得 
To ctp 
gss a i = a (126) 
根据 雷达 方程 , 此 时 对 应 的 回 波峰 值 功率 为 
P,G? Xo? A 
P= ve 


式 中 , 4 为 电磁 波 的 照射 面积 . 根据 Walsh 公式 (41, 
脉冲 有 限 足 迹 直 径 和 照射 面积 分 别 为 


D=2)/Hoer(1+ 5"), (18) 
4=7(2) = rHocr(1 + 2), (19) 


上 式 中 考虑 了 火星 的 曲率 ，Rm 为 火星 半径 . 将 式 
(16) 和 (19) 代入 式 (17) 可 得 


P,G? *er|Po|? Ho ctp 
r=  g(4n)?o2 HS (1 is Ee) 2 一 ) 
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这 就 是 粗糙 表面 的 高 度 计 雷达 方程 ， 该 方程 可 用 
于 HF 高 度 计 地 表 回 波 功率 的 提取 . 


3 火星 次 表层 介质 分 层 介 电 特 性 仿真 


一 般 地 , 火星 表面 和 次 表层 介质 都 可 以 用 空气 、 
水 ( 冰 )、 土壤 等 的 混合 物 建 模 . 在 高 于 100kHz 的 频 
Bi, 该 混合 介质 的 介 电 常数 可 用 下 式 表 示 M: 


Em = Ells) Etsel, (21) 


AF, em 为 混合 物 的 介 电 常 数 ; sa 为 空气 介 电 常 数 ， 
非常 近似 1; ew 为 水 的 介 电 常数 ; ss 是 固体 土壤 的 介 
EAR AF o 为 孔 性 ; 而 水 的 各 种 形态 (OK. WS 
纯净 水 和 盐水 等 ) 在 孔 性 中 所 占 比例 称 为 饱和 度 , 记 
H s ( 见 图 3). 根据 Debye 方程 66, 有 


Ew 二 E06 + = 一 ee 一 
1 — jw7, 
E 十 (Est = £e) j WT; (Est F E00) 
Tea 1+ (wr)? ¢ (22) 


火星 的 表层 和 次 表层 介质 中 富 含 碎 石 、 瓦 砾 等 . 
这 些 碎片 都 是 散射 体 , 使 得 介质 内 部 有 更 复杂 的 介 电 
结构 . 如 果 对 这 些 雁 片 边界 的 散射 行为 逐一 考查 , 整 
个 散射 行为 将 过 于 复杂 ， 以 至 于 无 法 计算 只 要 碎 
片 (散射 体 ) 在 介质 中 均匀 分 布 , 即 可 以 把 这 些 散 射 
体 和 背景 介质 等 效 成 单一 介质 . 

Kong 推导 了 球形 散射 体 的 等 效 介 电 常数 51: 
1 十 2 人 Ke 
1— foke 


oliz 


ĉe =e» | + jf (27)*- 


(23) 


2 (1-A) 
G+ 这 | 
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其 中 ， 
Ed — Eb 
ke = 24 
E gd + 2ep (24) 


Eb 为 背景 介 电 常 数 , sa 为 散射 体 (碎片 ) 的 介 电 常数 ， 
Fy 为 雄 片 的 体 占 比 , a 为 散射 体 半径 ，》 为 电磁 波 波 
K- 式 (23) 虽然 是 针对 球形 散射 体 推导 的 , 但 对 于 大 
体 光 滑 的 一 般 形 状 碎片 也 是 适用 的 . 火星 表面 的 碎片 
基本 是 随机 分 布 的 , 因此 可 用 式 (23) 对 含 碎 片 的 介 
质 介 电 常 数 进行 建 模 . 

根据 火星 特点 , 仿真 模拟 了 两 种 情形 下 的 介质 分 
层 模 型 , 并 计算 了 每 层 的 介 电 常 数 实 部 和 损耗 正切 . 

情形 1 浅 层 渗流 .其 160m 以 下 的 水 相 是 盐 
水 , 可 以 看 到 , 盐水 成 分 的 损耗 正切 很 大 , 在 此 情况 
下 电磁 波 的 穿 透 深度 受到 很 大 局 限 . 

情形 2 ARR. BHA AZE 50 
米 以 下 有 一 层 很 深 的 冰 透 镜 体 , 表现 为 纯 冰 . 其 损耗 
正切 比分 凝 冰 更 小 (0.0012), 电磁 波 在 其 中 的 穿 透 深 
度 最 深 ， 冰 透镜 体 和 分 凝 冰 都 是 火星 探测 中 很 重要 
的 固态 水 探测 要 素 . 

仿真 结果 列 于 表 2 和 表 3. 表 中 给 出 了 各 层 的 孔 
VE. 饱和 度 、 水 相 (水 相 为 0 代表 空气 , 1 代表 冰 , 2 代 
表 液 态 水 , 3 代表 盐水 )、 介质 材料 等 的 参数 作为 输入 
参数 . 输出 参数 为 介 电 常 数 实 部 和 损耗 正切 . 在 仿真 
中 , 假定 大 部 分 介质 的 铁 氧 体 含量 为 10%, 直径 0.1m 
的 碎片 体 占 比 为 5%. 


4 回 波 仿真 结果 


为 了 验证 火星 次 表层 分 层 介质 模型 和 反 演 方法 ， 
对 以 上 两 种 情形 的 回 波 进 行 了 仿真 . HF 高 度 计 对 各 
层 的 分 辨 率 如 图 4 所 示 . 由 图 4 中 可 见 , 各 层 的 分 辨 


aam 
| | air(é,=1) 


| water/ice, saturation 


BAX (s, % volume) 

E soil, porosity 

se (Q, % volume) 

Y iron (conductive material), 
WA concentration (C) 


图 3 火星 表层 内 含 空气 、 水 OK), 土壤 的 介质 模型 


Fig.3 Mars media model containing air, water (ice) and soil 
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(a) 


15 


resolution/m 


P 9 5 ? 
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layer index layer index 
图 4 各 层 垂直 分 辩 率 
Fig.4 Vertical resolutions of each layer 
表 2 情形 1 介 电 常数 仿真 结果 (20 层 ) 
Table 2 Dielectric constant simulation results of Case 1 (20 layers) 

层 数 ” 距 表面 高 度 /m 和 孔 性 /(%) ”饱和 度 /(%) KE 。 介质 材料 介 电 常数 实 部 ”损耗 正切 
1 1 50 0 0 风 沉 积 层 2.9922 一 0.0053 
2 3 15 0 0 硬 沉积 层 5.9524 一 0.0086 
3 5 50 0 0 沉积 玄武 岩 2.9922 一 0.0053 
4 10 5 0 0 致密 玄武 岩 7.2482 一 0.0095 
5 100 10 0 0 分 层 玄 武 岩 6.5680 一 0.0091 
6 110 50 100 1 风 沉 积 层 5.1839 一 0.0058 
7 150 10 50 1 分 层 玄武 岩 6.9401 —0.0091 
8 152 20 0 0 冲积 沉积 层 5.3951 一 0.0081 
9 160 50 0 0 火山 灰 2.9922 一 0.0053 
10 200 10 100 3 分 层 玄武 兰 12.9450 一 0.1544 
11 220 20 100 3 ”火山 口 溅 射 物 20.5938 一 0.3089 
12 250 10 20 3 分 层 玄武 岩 7.5303 一 0.0376 
13 255 50 20 3 风 沉 积 层 5.8696 一 0.1484 
14 350 10 20 3 BERS 7.5303 一 0.0376 
15 355 20 20 3 冲积 沉积 层 7.0831 一 0.0651 
16 500 10 20 3 分 层 玄武 岩 7.5303 一 0.0376 
17 750 10 20 3 分 层 玄武 岩 7.5303 一 0.0376 
18 760 50 20 3 火山 灰 5.8696 一 0.1484 
19 900 10 20 3 DEERE 7.5303 —0.0376 
20 1000 10 20 3 分 层 玄 武 岩 7.5303 一 0.0376 


率 位 于 5~10m. 各 层 的 衰减 如 图 5 所 示 , 可 以 看 到 ， 
情形 2 中 第 6 层 ( 冰 透 镜 体 ) 的 衰减 远 比 其 他 层 低 ( 约 
为 0.01dB.m-1), 因此 高 度 计 电磁 波 在 冰 层 中 可 穿 透 


根据 式 (2) 可 求 出 每 两 层 之 间 反 射 回 波 对 应 的 
EIR, 对 式 (6) 和 式 (7) 进行 推广 , 可 得 各 层 之 间 界 
面 反射 回 波 对 应 的 回 波 功率 . 将 这 些 延迟 与 功率 对 应 
起 来 可 得 到 回 波 波 形 , 如 图 6 所 示 . 纵 坐 标的 功率 值 


N 
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表 3 情形 2 介 电 常数 仿真 结果 (8 层 ) 
Table 3 Dielectric constant simulation results of Case 2 (8 layers) 
层 数 。” 距 表面 高 度 /m ”了 和 孔 性 /(%) ”饱和 度 /(%) 水 相 。 介质 材料 介 电 常数 实 部 ”损耗 正切 
1 1 50 0 0 风 沉 积 层 2.9922 —0.0053 
2 3 15 0 0 便 沉 积 层 5.9524 —0.0086 
3 5 50 100 1 沉积 玄武 岩 5.1839 一 0.0058 
4 50 10 100 1 ”致密 玄武 岩 7.3336 一 0.0092 
5 55 50 100 1 RZA 5.1839 一 0.0058 
6 205 100 100 1 冰 透 镜 体 3.2000 一 0.0012 
7 225 20 100 1 冲积 沉积 6.7234 一 0.0083 
8 1000 10 100 1 4 anne 7.3336 —0.0092 
0.50 r 0.09 -一 -一 一 — Tee 
0.45 “全 l (a) 0.08 Case 2 (b) | 
一 0.40 ZS WO 
POA E 0.06 
ə 0.30 2 0.05 
Ate E 0.04 J 
§ 0.20 3 
he 3 0.03 
3 910 = 0.02 
0.05 0.01 g 
i ef f J 
5 2 25 1 2 3 4 5 6 7 8 9 
ka toe layer index 
图 5 各 层 的 衰减 
Fig.5 Attenuations of each layer 
0 OO 0 r — 
Case | (a) > 
-50 i i -50 Case 2 HE 
一 100 4 —100 
a -150 | 各 -150 
5 -200 | 5 -200 
2. -250 J 2 250 1 
= = 
2 -300 4 5 -300 ] 
-350 J -350 J 
-4001 a as a -400 eet | 3 
02 04 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 («10*) 0 0.2 04 0.6 0.8 10 12 1.4 1.6 1.8 2.0 (x10^ 
delay/ns delay/ns 
图 6 各 层 回 波 波形 (冲击 脉冲 ) 
Fig.6 Echo waveform of each layer (impulse) 
已 转换 为 dBm， 横 坐标 利用 系统 的 压缩 脉冲 宽 图 6 的 仿真 中 将 高 度 计 脉冲 看 为 冲击 脉冲 , 没有 


BE (1/10MHz=100ns) 进行 了 归 一 化 , 对 应 于 回 波 ”考虑 


的 时 间 分 辩 单 元 . 


雷达 系统 点 日 


标 响应 的 景 


CMe. 无 限 光滑 平面 的 


反射 行为 可 看 做 点 目标 , 而 在 高 度 计 系 统 中 使 用 了 


RUS F: 基于 分 层 介质 模型 的 HF 雷达 


echo power/dBm 


—3500 4 
0 0.2 04 0.60.8 101.2 14 1.6 1.8 2.0(«104) 
delay/ns 
图 7 各 
Fig. 7 


全 去 斜 脉冲 压缩 技术 , 其 点 目标 响应 对 应 于 since K 
BO, 通常 为 了 数学 表达 和 推导 的 简便 ， 可 将 其 表 
示 为 高 斯 函数 形式 , 即 


2 


S,(t) = Aexp ( — za) (25) 
st, 4 为 归 一 化 幅度 ; op = k7, 7 为 系统 的 压缩 脉 
冲 宽度 (带宽 的 倒数 , 100 ns); k 为 点 目标 响应 等 效 高 
斯 函数 的 比例 因子 , 根据 Jason 高 度 计 团 队 的 研究 结 
果 , 取 k= 0.513128), 

因此 , 高 度 计 接 收 到 的 总 的 回 波 功率 

=R 
202 


P(t) => Bap] |; 
i=0 
式 中 , n 为 介质 的 层 数 , P., 为 第 i 层 介质 底部 的 反 
射 功率 , 为 第 i 层 介质 底部 回 波 对 应 的 延迟 . 
根据 式 (26) 得 到 的 高 度 计 回 波 波形 如 图 7 所 
示 ， 从 图 6 和 图 7 可 以 看 到 , 在 情形 1 中 , 由 于 次 
表层 200m 左右 的 介质 介 电 常数 过 大 ( 实 部 大 于 10, 
虚 部 接近 2), 电磁 波 急剧 衰减 , 更 深 的 次 表层 难以 探 
W. 而 在 情形 2 中 , 由 于 系统 分 辩 率 的 限制 , 1~3 层 
和 4~5 层 实际 上 表现 为 1 个 峰 , 冰 层 位 于 第 10~25 
个 回 波 单元 之 间 ， 冰 层 以 下 的 反射 功率 明显 降低 很 
多 . 这 是 因为 冰 层 以 下 的 反射 系数 明显 降低 , LER 
增高 ， 因 此 通过 回 波 下 降 的 幅度 和 持续 的 延迟 可 识 
别 冰 层 的 分 布 
以 上 讨论 的 是 理想 光滑 分 层 界面 的 回 波 , 当 表 面 
或 各 层 介质 之 间 的 界面 为 粗糙 时 , 就 要 考虑 如 式 (10) 
所 示 的 粗糙 度 等 参数 , 实际 的 回 波 功率 会 降低 . 但 这 


(26) 
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echo power/dBm 


层 回 波 波形 (实际 点 


Echo waveform of each layer (actual PTR) 
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—3500 m 
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delay/ns 


目标 响应 ) 


里 仿真 的 是 单 脉冲 回 波 功率 , 对 于 多 脉冲 回 波切 率 还 
可 以 进行 累积 . 对 于 传统 模式 , 回 波 功 率 的 累积 是 非 


相关 的 , 回 波 功率 要 乘 以 累积 脉冲 数 的 平方 根 ; 对 于 
合成 孔径 模式 , 累积 是 相关 的 , 回 波切 率 要 乘 以 累积 
脉冲 数 . 

5 结论 


介绍 了 HF 雷达 高 度 计 系统 的 原理 和 参数 设计 
推导 了 高 度 计 电磁 脉冲 和 多 层 光 滑 介质 之 间 的 相互 
作用 , 并 提取 其 回 波 功 率 随 时 间 变 化 的 模型 , 同时 考 
虑 了 表明 粗糙 度 对 表面 回 波 功率 的 影响 ， 选 用 两 种 
典型 的 火星 分 层 介质 模型 , 并 对 其 各 层 的 介 电 常数 特 
性 和 高 度 计 回 波 时 间 延 迟 、 功 率 衰 减 和 波形 形状 等 
进行 了 仿真 . 仿真 结果 表明 , 采用 HF 雷达 高 度 计 系 
统 可 以 实现 对 火星 次 表层 介 电 特 性 垂直 廓 线 的 反 演 ， 
从 而 实现 对 火星 中 的 水 (包括 固态 水 和 液态 水 ) 的 识 
别 和 探测 . HF 高 度 计 对 于 火星 次 表层 的 探测 具有 有 
重要 意义 , 在 今后 中 国 自 主 火星 探测 中 可 考虑 选用 该 
ASB. 
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